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Abstract 

In this paper, we study the Lie point symmetry group of a system describing an ideal plastic 
plane flow in two dimensions in order to find analytical solutions. The infinitesimal generators that 
span the Lie algebra for this system are obtained. We completely classify the subalgebras of up to 
codimension two in conjugacy classes under the action of the symmetry group. Based on invariant 
forms, we use Ansatzes to compute symmetry reductions in such a way that the obtained solutions 
cover simultaneously many invariant and partially invariant solutions. We calculate solutions of the 
algebraic, trigonometric, inverse trigonometric and elliptic type. Some solutions depending on one or 
two arbitrary functions of one variable have also been found. In some cases, the shape of a potentially 
feasible extrusion die corresponding to the solution is deduced. These tools could be used to thin, 
curve, undulate or shape a ring in an ideal plastic material. 

Running Title: Symmetry group analysis of an ideal plastic flow 
PACS numbers: Primary 62.20.fq; Secondary 02.30.Jr 

Keywords: symmetry group of partial differential equations, symmetry reduction, invari- 
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I Introduction 

In this paper, we investigate the plane flow of ideal plastic material^^ modelled by the 
hyperbolic system of four partial differential equations (PDE) in g = 4 dependent variables 
0", 6, u, V and p = 2 independent variables x and 

(a) - 2k {6^ cos 26 + 9y sin 29) = 0, 

(b) ay ~ 2k {e^sm2e -Oy cos 20) = 0, 

(c) {uy + Vx) sin 29 + {ux — Vy) cos 26 = 0, 

(d) Ux + Vy = 



'1 



where cr^ = da/dx, etc. The expressions ([Ha), ([Hb), are the equilibrium equations for 
the plane problem. In other words they are the Cauchy differential equations of motion 
in a continuous medium where we consider that the sought quantities do not depend 
on z. These two equations involve the dependent variables a and 9 that deflne the 
stress tensor; a is the mean pressure and 9 is the angle relative to the x axis in the 
counterclockwise direction minus 7r/4. The equation ([He) corresponds, in the plane case. 
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to the Saint-Venant-Von Mises plasticity theory equations, where u and v are respectively 
the velocities in the x axis and y axis directions. Moreover, we assume incompressibility 
of the material and hence the velocity vector is divergenceless. This explains the presence 
of the equation (fUd) in the considered system. The positive-definite constant k is named 
the volumetric compression coefficient and is related to the Poisson coefficient and the 
Young modulus v by the formula k = {1 — 2v)E~^. 

In order to calculate new solutions of the system consisting of (IHa) and (fUb), S.I. 
Senashov et al^^ acted with transformations of the symmetry group of this system on 
known solutions of some boundary value problems, i.e. the Nadai solutionPfor a circular 
cavity under normal stress and shear and the Prandtl solutioi]^ for a bloc compressed 
between two plates. In addition, Czy^ found simple and double wave solutions for the 
system ([1]) using the method of characteristics. However, as it is often the case with 
this method, his solutions rely on numerical integration for obtaining the velocities u and 
V. To our knowledge, no systematic Lie group analysis based on a complete subalgebra 
classification in conjugacy classes under the action of the symmetry group G of the system 
([1]) has been done before. 

The goal of this paper is to systematically investigate the system ([1]) from the per- 
spective of the Lie group of point symmetries G in order to obtain analytical solutions. 
That is, we obtain in a systematic way all invariant and partially invariant (of struc- 
ture defect 5 = 1 in the sense defined by Ovsiannikov^) solutions under the action of G 
which are non-equivalent. Invariant solutions are said non-equivalent if they cannot be 
obtained one from another by a transformation of G (the solutions are not in the same 
orbit). In practice, we apply a procedure developed by J. Patera et al^^ that consists 
of classifying the subalgebras of C associated with G into conjugacy classes under the 
action of G. Two subalgebras Ci d C and C[ d C are conjugate if GCiG = C'^. For each 
conjugacy class, we choose a representative subalgebra, find its invariants and use them 
to reduce the initial system ([Tj) to a system in terms of the invariants which involve fewer 
variables. According to the approach proposed by Kruskal and Clarkson,^^ which is part 
of the more general framework of conditional symmetries, we propose in this paper some 
Ansatzes which allow us to cover simultaneously many invariant and partially invariant 
solutions (PIS). These more general solutions reduce to invariant and partially invariant 
ones for appropriate parameter values. We illustrate these theoretical considerations with 
many classes of solutions. A more exhaustive collection of such solutions is provided in^^ 
Thereafter, we draw for some solutions the shape of the corresponding extrusion die. The 
applied method relies on the fact that the walls of the tools must coincide with the flow 
lines describe by the velocities u and v of the solutions of the problem. In application, 
it is convenient to feed in material the extrusion die rectilinearly at constant speed. So, 
the tools illustrated in this paper were drawn considering this kind of feeding. Based on 
mass conservation and on the incompressibility of the materials, we easily deduce that 
the curve deflning the limit of the plasticity region for constant feeding speed must obey 
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the ordinary differential equation (ODE) 

dy ^ Vo-v{x,y) 
dx Uq — u{x, y) ' 

where Uq, Vq are components of the feeding velocity of the die (or extraction velocity at 
the output of the die) respectively along the x-axis and y-axis. One should note that the 
conditions are reduced to those required on the limits of the plasticity region in the 
paper of Czy2Pwhen Vq = and that the curves defining the limits coincide with slip lines 
(characteristics), that it when we require dy/dx = ta.n6{x,y) or dy/dx = — cot6{x,y). 
Thus the condition ([2]) can be viewed as a relaxation of the boundary conditions given in 
the work of Czyz.^ The reason we can use these relaxed conditions is that we choose the 
walls of the tool to coincide with the flow lines for a given solution rather than require 
the flow of material to be parallel to the walls. Using these relaxed conditions, we can 
choose (in some limits) the feeding speed and direction for a tool and this determines the 
limits of the plasticity region. 

The paper is organized as follows. In section |Tl] we give the infinitesimal generators 
spanning the Lie algebra of symmetries C for the system ([1]) and the discrete transfor- 
mations leaving it invariant. A brief discussion on the classification of subalgebras of C 
in conjugacy classes follows. Section IIIII is concerned with symmetry reduction. It de- 
scribes how the symmetry reduction method (SRM) has been applied to the system ([T]) 
and the method for finding partially invariant solutions. More precisely, we give several 
results obtained from Ansatzes so that each presented solution includes many invariant 
and partially invariant solutions corresponding to appropriate choices of the parameters. 
We conclude this paper with a discussion on the obtained results and some incoming 
results. 



II Symmetry algebra and classification of its subalgebras 

In this section we study the symmetries of the system ([1]). Following the standard algo- 
rithm,^^ the Lie symmetry algebra of the system has been determined. It is spanned by 
the eight infinitesimal generators 

Di = xd^ + ydy, D2 = udu + vd^, B = -yd^ + xd^, 
Pi = dx, P2 = dy, P3 = dfj, P4 = du, P5 = dy, 

where we use the notation dx = d/dx, etc. The generators Di and D2 generate dilations 
respectively in the space of independent variables {x, y} and the space of dependent 
variables {u, v}. Moreover, B is associated with a kind of boost and the Pj, z = 1, . . . , 5 
generate translations. The commutation relations for the generators are shown in table 
m One should note that the system ([T]) is invariant under the discrete transformations: 

-Ri : X —X, y 1— i- —y, a a, 9 9, u ^ u, v v; 
R2 ■ X X, y y-^ y, a ^ a, 9^9, —u, v (-)■ —v. 
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Table I: Commutation relations for the algebra £. 
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These reflections induce the automorphisms of the Lie algebra C: 



7^l 






D2 ^ D2, 


B ^ -B, 


Pi ^ -Pi 




P2 


^-P2, 


P3 ^ P3, 


Pa ^ P4, 


A ^ P5; 


7^2 






D2 ^ D2, 


B ^ -B, 


Pi ^ Pi, 




P2 


^P2, 


Pz ^ Pz, 


P4 ^ -P4, 


P5 ^ -A- 



(5) 



Since we seek solutions that are invariant and partially invariant of structure defect 5 = 1, 
we only have to classify the subalgebras of codimension 1 and 2. We have used the 
following factorization of the Lie algebra C: 

C = {{{D^, D2} > {B}} > {Pi, P2, P4, A}} © {P3} , (6) 

where > denotes the semi-direct sum and © the direct sum of Lie algebras. The Lie 
algebra C is the direct sum of the center {P3} with the subalgebra {{{Di, D2} > {B}} > 
{Pi, P2, P4, P5}} which contains an abelian ideal {Pi, P2, P4, P5}. Applying the method,'^^ 
we proceed to classify all subalgebras of C in conjugacy classes under the action of the 
automorphisms generated by G and the discrete transformations (jlj). In practice, we can 
classify the subalgebras under the automorphisms generated by G and decrease the range 
of the parameters that appear in the representative subalgebra of a class using the Lie 
algebra automorphisms (|5]). The classification results are shown in table [IV] for subal- 
gebras of codimension 1 and in tables IIIII and IIVI for subalgebras of codimension 2. For 
each subalgebra in these tables, a complete set of invariants is given. Table IIIII contains 
the codimension 2 subalgebras which admit a symmetry variable, while table IIVI lists the 
codimension 2 subalgebras that have no symmetry variable. Despite the absence of a 
symmetry variable, the subalgebras in table [IV] lead to PIS. 

Ill Symmetry reductions and solutions of the reduced systems. 

In this section we use the symmetry reduction method, as presented in,'^ to compute 
invariant solutions under the action of subgroups of the symmetry group G of the initial 
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system ([T]). Following the usual reduction procedure!^ we consider a subgroup Gdj C G 
associated with a subalgebra Cdj C £ of dimension d. Then the subgroup Gdj admits 
k = p + q — d functionally independent invariants / = {Ii{x,u), . . . ,Ik{x,u)), where x 
denotes the for independent variables and u the dependent variables. We have to consider 
three different possibilities. 

(i) We first make the following hypotheses: 

• q < d < p + q, 

• the rank of the Jacobian matrix {dl/du) is q, i.e. Ta.nk{dl /du) = q, 

• the complete set of invariants under the action of Gdj takes the form 

{C,i{x), . . . ,^k-q{x), Ii{x,u), . . . , Iq{x,u)}, whcrc denotes the symmetry vari- 
ables. 

With these conditions, the invariant solutions under Gdj are in the form 

= U,{x, ^ = 1, . . . , g, F(0 = . . . , FqiO), (7) 

where U = {Ui, . . . , Uq) is the solution u = {ui, . . . , Uq) of the relations 

F,{0 = h{x,u), ^ = {Ux),...,^k-,). (8) 

The Fi functions that depend on the k — q variables ^j, satisfy the differential equation 
system that results from the introduction of ([7]) into the system ([1]). To calculate the 
invariant solutions of the system ([1]), we must consider subalgebras of codimension 
d = 1. Indeed, the system ([1]) is expressed in term of p = 2 independent variables, so 
we cannot reduce it by more than one variable. In other words, we must have k = 5 
to apply the SRM to the system ([1]). This is possible only ii d = 1. 

(ii) Consider now the situation where Ta.Bk{dI /du) = q' < q and {C,i{x),..., ^k-q'i^), 
Ii{x,u), . . . ,Iq/{x,u)} is a complete set of invariants of Gdj- In this case, it is only 
possible to solve the relations 

F,(0 = /,(x,m), j = l,...,g', e = (6,---,6-a')- (9) 

for q' of the dependent variables u. Without loss of generality, we suppose that 
we solve the relations for the first q' dependent variables {ui, . . . ,Uq>). Then, 
introducing 

= U,{x, ^ = 1, . . . , g', F(0 = . . . , FqiO), (10) 

into the system ([1]), we find a system of differential equations for the functions Fi{^) 
in which we have the remaining Uqi^i{x), . . . ,Uq{x), which depend on the original 
independent variables x. We must add the compatibility conditions on the mixed 
derivatives of Wg'+i, . . . ,Uq. The obtained solutions are partially invariant. In this 
paper, we are interested in PIS of structure defect 6 = 1 which are obtained from 
codimension 2 subalgebras. Subalgebras which satisfy the condition rank(9//9'u) = 
3 < g are listed in table IIIII 
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(iii) In the case where the subgroup Gdj, corresponding to a codimension 2 subalgebra 
Ldj, has a complete set of invariants satisfying the condition rank ((9//(9u) = 4, 
there always exists, for the considered system ([1]), an invariant in terms of x, y and a 
only. We denote this invariant r. Consequently, the set of invariants takes the form 
{r(x, y, a), li{0, u, v), hiO, u, v), I-s^O, u, v)} and we can solve for 9, u, v the relations 

F,ir)=I,ix,u), 1 = 1,..., 3, (11) 

and introduce the result in the system ([T]). To justify this particular choice for the 
invariant r associated to Cdj, one should note that the equations ([Ua), ([Hb), do not 
involve the quantities u,v and that the quantity 6 is an invariant of the whole Lie 
algebra C. This choice leaves the equations ([Ha), ([Hb), uncoupled to the equations 

dHc), dHd). 



In the case of (i) and (ii), we can always choose the invariant Ii = 9. Therefore, from 
the relations ([9]) we have 

e{x,y) = J{a^,y)), (12) 

where we write Fi as J for simplification. We will keep this notation for the subsequent 
sections of the paper. Introducing 9 given by (fT2|) in the equations ([Ha), ([Hb), and using 
the compatibility conditions on the mixed derivatives of a relative to x, y, we obtain 
an ODE for the function J. If we find the solution for J, the compatibility condition is 
satisfied and then we can integrate by quadrature to find the solution for a. This solution 
is not necessarily an invariant one, but it includes invariant solutions for appropriate 
choices of parameters that appear in a. 



Concerning the case (iii), since there is no symmetry variable, we propose the solution 
for 9 in the form 

9{x,y) = J{T{x,y,a{x,y))). (13) 

This allows one to solve the equations ([Ha), ([Hb), for 9, a, without considering the 
equations ([He) and ([Hd), because this leaves the equations uncoupled. None of the 
codimension 2 subalgebras of table jlVj admit a symmetry variable, but they possess an 
invariant in the form r = hi{x,y)a + h2{x,y). Thus, we can express a as 

a = L^J^ (14) 
hi{x,y) 

and eliminate the quantities a from the equations ([Ha), ([Hb), and then use the compat- 
ibility condition on mixed derivatives relative to x and y of the function r to obtain an 
ODE that must be verified by J. If we find the solution for J, we can find r by quadrature 
and consequently we find a through the relation (I14p . 
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Below we present, for each distinct form of the symmetry variable ^ or of the invariant 
r, the most general solution for 9 and a of the equations ([Ha), (IHb), where we suppose 
a solution for 9 given by (I12p or (fT5]) depending on the considered case. Then for each 
of these solutions, we consider Ansatzes on the form of the velocities u and v to compute 
solutions that cover all invariants solutions corresponding to a given form of the symmetry 
variable ^ or of the invariant r. In the table |Vl the subalgebras £j j are classified by 
symmetry variable and according to the form of the suggested Ansatz on the form of the 
solutions for u and v. In general, the solutions are not necessarily invariant under the 
action of the symmetry group G but they are so for an appropriate choice of integration 
constants. We illustrate with some examples the suggested method. 



A Solution for 9 in the form of a propagation wave. 

The goal of this section is to construct a solution for the angle 9 in the form of a prop- 
agation wave and to get the corresponding solution for the pressure a and then, in the 
subsections [H [21 we get the velocities u and v from different assumptions made on the 
form of their solutions, i.e. we make additive separation and a multiplicative separation. 



We look for a solution of the system ([Ha), ([Hb), such that 9{x, y) has the form 

9{x,y) = J{ax,y)), (15) 

with the function J to determine and where 

^{x,y) = aix + a2y, ai,a2GE. (16) 

This type of solution includes, for appropriate values of the parameters ai, a2, all invariant 
and partially invariant solutions (for the quantities 9, cr) corresponding to subalgebras Lij, 
j = 1, . . . , 4 and L2^k, k = 1,4, 5, 6, 7, 8, 10, 11, (that is the subalgebras listed at line no. 1 
in table [Vl). Introducing (llSp . ([TB]) in ([Ha), ([Hb), we get that the following system must 
be verified 

(a) a,{x, y) = 2k [ai cos{2J{^{x, y))) + as sm{2J{^{x, y)))] J'i^{x, y)), 

(b) ay{x, y) = 2k [ai sin(2 J(e(a;, y))) - oo^{2J{^{x, y)))] J\i{x, y)). ^ ^ 

The compatibility condition on the mixed derivatives for cr relative to x and y provides 
the next ODE for J in term of ^ 

[-A sin(2 J(0) + /i cos(2 J(0)] J'\i) + [-2^^ sin(2 J(0) - 2A cos(2 J(0)] J'{if = 0, (18) 
where Ci,C2 G R are integration constants. The solution of f ll8p is 



J(0 = i^vr + 1/2 arctan ( f^i^^^ + + ' i^^^ + , , (ig) 



where s denote the sign of the expression A(ci^ + C2) — + jj'^ — (ci^ + 02)^ and 

\ = a\ — a\, IX = —2aia2- The solution for 9{x,y) is provided by ([T5l) with J defined by 



7 



T9|) and ^{x,y) by f|T6|) . Thereafter, we substitute the solution (fT9l) for J in the system 
T7|) and we solve for a by quadrature. We find the solution 



, , cik(a2X-aiy) J (aj + aj) - {ci{aix + a2y) + 
<y{x,y) = — — 2 — 2 + ^3, C3 G R. (20) 

We are interested in real solutions, so 9 and a are defined over the domain 

Vt = {{x,y) : (al + al) - (ci(aiX + aay) + Ca)^ > O} . (21) 

Thus, we have that the quantities 6 and a of all invariant solutions for the above sub- 
algebras are covered by 9{x,y) = J{C,{x,y)) with J defined by (1191) . ^{x,y) by (ITB]) and 
(T{x,y) by m- 

1 Additive separation for the velocities u and v. 

Corresponding to the solution for 6 and a given respectively by (fT5|) . (fT6l) . (fT9|) and by 
(120|) . we seek for a solution u and v of the additive separated form 

u{x,y) = f{x,y) + F{^{x,y)), v{x,y) = g{x,y) + G{^{x,y)). (22) 

The first step is to classify the admissible forms for the functions / and g. To do this, we 
substitute (l22l) in the system ([Uc), ([Hd), and we obtain the new system 

(a) [/, + a2F'(e(x, y)) + + a,G'{a^, y))] sin(2 J(e(x, y))) 

+ [f. + aiF'i^x, y)) -gy- a2G"(e(x, y))] cos(2 J(e(x, y))) = 0, (23) 
(6) + aiF'(e(x, y)) + gy + a2G'{i{x, y)) = 0, 

where F', G' represent the derivative of F, G relative to ^ respectively. We will keep this 
notation in all of the following sections. We next introduce the linear differential operator 

= -a2-Pi + aiP2 = -a2dx + aidy, (24) 

which annihilates any function in term of the quantity C,{x, y) = aix + a2y. Now, applying 
it to the system (l23l) we find the following conditions that constrain / and g: 



(a) {L^ify + g.,)) sin(2 J(e(x, y))) + 2{L^Q cos(2J(e(a;, y))) = 0, 

(6) L^{U + gy)=Q. ^ ^ 

Then, taking into account that Lg, dx and dy commute with each other and using the 
notation 

f = L J, = L^g, (26) 
we can rewrite the conditions (1251) in the new form 



[a] (f, + g,) sin(2 J) + 2f, cos(2 J) = 0, 

{h) fx + S, = 0. ^ ' 
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We have to consider two cases separately, that is when f and g verify the equation (l271 a) 
by requiring that the coefficients of the trigonometric functions vanish 

U = fy + Q. = 0, (28) 
or by imposing that f and g obey the relation 

fy + 0x = -2cot(2J(0)fx. (29) 

(i) If we suppose that f|27[ a) is satisfied by the conditions f l2S]) . then we have the following 
relations 

(a) LJ, = 0, 

(6) L^{fy + g,) = 0. ^ ' 

If ai 7^ 0, 02 7^ 0, the solution for / and g of (130|) takes the form 

fix,y) = CiiO + ^iy^ 9ix,y) = (2(0 + ^2X, e R, 2 = 1,2,3. (31) 

where the functions Q are arbitrary functions of ^. If ai = 0, then we can set 02 = 1 
without loss of generality. In this case, the functions / and g take the form 

f{x,y) = uix + Ciiy), g{x,y) = U2X + C2iy), w^gR, z = l,2. (32) 

where the function Q are arbitrary functions of y. Similarly, we can consider the case 
Oi = 1, 02 = 0, which leads to the solutions 

fix,y) = Ciix) +uiy, g{x,y) = C2ix) +U2y, w^gR, 2 = 1,2, (33) 



where the functions Q are arbitrary functions of x. Under the hypothesis fl28|) . and 
according to appropriate values of a^, the form for the functions / and g given by fl3T]) . 
(152]) or fl33p . are the most general such that the Ansatz (122]) reduces the system ([Uc), 
([TJd), for M, f to a system of two ODE for F and G in term of ^ given by 

^ Qi(/x -^;/) + Q2(/y + ^x)sin(2J(0) + 2a2/xCos(2J(0) 

(a2-a2)sin(2J(0) + 2aia2Cos(2J(0) ' (34) 

G'{0 = -ci2HciiF'{0 + f. + 9y)- 

We distinguish two different cases for the solutions F and G of ( 1M]) . that is when ai 7^ 
and 02 7^ or when 02 = 0. 

(a) If Oi 7^ and 02 7^ 0, we obtain 

^^^^ ^ _^ ^i(^i+^2)cos(2J(0) ^ a2C.isin(2J(0) ^ ^^^^ 

2ci Cl 



and 



^/^N > /.N fli 0^1 + c^2 COS 2 J aiUJi sm 2J W1+W2 . . 

<^(U = -(2(0 ^ + C5- (36) 

Cl Cl a2 
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Redefining the parameters ui, U2, the solutions u and v of ([He), (IHd), are then 



u{x,y) 
v{x,y) : 



, ai(c4 + C5) cos(2J(^(x,?/))) a2C4sin(2(J(^(x,?/)))) 
C4a; H ;:; \ h C6, 



2ci 

ai(c4 + C5) cos(2J(^)) aiC4 sin(2 J(^)) 



2a2Ci 



Cl 



Cl 

ai(c4 + C5)x 
02 



- c^y + C7, 
(37) 

where the q are integration constants. So, we obtain a solution of the system ([1]) by 
defining the angle 6 by (fT5|) . (fT9|) . the pressure a by (120|) and the velocities u and f by 
(137|) . with ^{x,y) = aix + 02?/. As can be seen in the illustration, we plotted in figure 
[T] the shape of an extrusion die corresponding to the solution (137|) for the velocities u 
and V. The vector field defined by these velocities is traced inside the tool. The mean 
pressure in the extrusion die, which is given by (120|) . is represented by shades of grey, 
from pale grey for lower pressure to dark grey for higher pressure. The parameters have 
been set to k = 0.0027, ai = 1, 02 = 1, Ci = 0,5, C2 = 0, C3 = 0, C4 = 0, C5 = — 1, 
Ce = 0, C7 = and the feeding velocities modulus Uq, Ui, are both equal to ^/2, while the 
extraction velocities modulus Vq, Vi are both equal to a/3. The curves Ci, C2 determine 
the limit of the plasticity region at the mouth of the die while the curves C3, C4, have 
the same meaning at the end of the die. It's a double feeding tool which has a symmetric 
configuration under the reflection x 1— )■ —x except near the plasticity limits defined by 
C3 and C4. Other feeding and extraction settings are possible, that is the angle and the 
velocities modulus, but no setting with completely vertical extraction speeds were found. 




Figure 1: Extrusion die corresponding to the solution (PO)) . (|37p . 
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(b) Now, if ai = and a2 = 1, then / and g are given by (132|) . In this case, we find 



Fiy) = -Ciiy) -^2y -^uji J cot 2 J (y) dy + ^^^^ 

Giy) = -uJiy - C2(y) + C5. 
By redefining Ui, U2, the solution may be written as 

u{x, y) = C4X — c^y — 2c4 J cot 2J{y)dy + C5, 
vi.x.y) = -c^y + C7, 



(39) 



where the q are real arbitrary constants. So, a solution of the system ([T]) consists of the 
angle 6 by (fTSj) . (fT9|) . the pressure a by (120|) and the velocities m and f by (139|) . with 
C,{x,y) = y. For this solution, we did not trace an extrusion die because of the similarity 
with the next case where a tool will be plotted. 

If Oi = 1 and 02 = 0, then the functions / and g are given by (133|) . We find in this case 
F{x) = -Ciix) - U2X + C4, 

f (40) 
G{x) = —uix — C2{x) + 2uj2 / cot 2J{x)dx + C5. 



Redefining wi, U21 the solution for u and f of ([T]c), ([TJd) is 

m(x, y) = -C^X + C4?/ + C6, 



(41) 

f (x, ?/) = — 040; + C5?/ + 2c5 / cot 2 J(a;)(ia; + C7, 



where the Cj are real arbitrary constants. The system ([T]) is solved by the angle 9 by fllSp . 
(fT9]) . the mean pressure a by (12(11) and the velocities u and v by (HTj) . with C,{x,y) = x. 
An extrusion die corresponding to the solution (1411) . with the mean pressure provided by 
(120|) . is shown in figure |2]for parameter values set to A; = 0.0027, Oi = 1, a2 = 0, Ci = 0, 5, 
C2 = 0, C3 = 0, C4 = — 1, C5 = 1, ce = 2, C7 = 0. The pressure is represented by the shades 
of grey ranging from pale to dark as the pressure increase. The curves Ci and C2 define 
the plasticity regions respectively at the mouth and at the end of the tool. The feeding 
velocity used is [/q = 1 and the extraction velocity is Vq = 1. For example, this die would 
be used to curve a plate or a rectangular rod. 

(ii) Now, we consider the case where f and g satisfy the condition (129|) . By the use of 
(p7|) . the compatibility condition on q leads to the hyperbolic equation 

-fyy-2 cot(2J(0)fx, + 4a2 csc(2J(0) J'(Ofx = 0, (42) 
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Figure 2: Extrusion die corresponding to the solution (PO)) . (j37p . 



over all the domain fl defined by (12T|) . So, we introduce the variables change {x,y) to 
{ip,ip) defined by 

ip = X -\ (ai cos(2 J) + 02 sin(2 J) — 2 J) , 

Cl 

i/j = x + — (ai cos(2 J) + 02 sin(2 J) + 2 J) , 

Cl 

which enables us to rewrite in the simplified form 



(43) 



U^iL„^_Weot,(2.) + ^i±|. (44) 

One should note that taking the difference of equations f l43p . we deduce 

J = S^Zll. (45, 

4a2 



We can find a particular solution of dHj) by assuming the solution in the separated form 

4o2 



f(V9, '?/') = P{K) + Q( J) where we use K = £1^±^. Under such a hypothesis, the solution 



of the equation (H3|) is 

f ^) = + ^) + a;2 cos f + Y ( ) + ^3, (46) 
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where the G R, i = 1, . . . , 3, are integration constants. By considering the change of 
variable ( 143|) . the system (1271) takes the form 

(^) sin (^) 

Sin (a^) cos 
from which we find the solution for g defined by 

V,) = -i^ cos f £1%^) + siu f £i<:p^l - + *) + -4. (48) 



2 V 2a2 / V / 2a2 



where G R, i = 1, . . . , 4. Finally, returning to the initial variables and considering the 
relations f l26p . we obtain / and g defined by 

ujici{2ai^{x,y) - a2ri{x,y))ri{x,y) [[ wiaA 

fix,y) = . ^ 2 I + + — cos 2J e 

4a2(ai + a^j V V 2a2 / 

+a;i sin(2J(^(x,?/))))r7(x,?/) + a;3?7(a;,?/) + Ci(^(a;,?/)), 

(49) 

where we introduced the notation 

^{x,y) = aix + a2y, r]{x,y) = -a2X + aiy 
and where 0(0 ^^^^ ^^^^^ arbitrary functions of one variable. 



It remains for us to integrate the equations fl23j) with / and g given by fl49l) . which form 
a system of two ODE for the unknown quantities F and G in term of ^. The solution for 
F is expressed by the quadrature 

F(0 = j [2a2K3 ((a2 - 1) sin(2 J) + 2ai cos(2J))]-' 

X uJi(^ - aici {ki sin(2J(^)) - 2a2 cos(2J(^))) ^ - K3 (2a2Ki sin^(2J(^)) 

+ (ai + 802)^2 cos(2 J(0) sin(2 J(0) - 2aia2 cos^(2J(0) 

+ a;2 ( - aiCiKi^ sin(2 J(0) - 2^3 cos(2 J(0) (k? sin(2 J(0) - 2a^ cos(2 J(0)) 
— 2a;3a2ft3 (ki sin(2 J(^)) — 2a2 cos(2 J(^))) + 2c(;4a2KiK3 sin(2 J(,^)) 



+ 2ft:3a2 ( («:2 sin(2 J(0) - 2ai cos(2 J(0)) ^ ' ^2 sin(2J(0)^) 



dC, + C4. 

(50) 
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It's the same for the solution for G, 
G(0 = J [2a,a2K3 {{a^ - 1) sin(2J(0) + 2ai cos(2 J(0))]"' 

oji (^aiCi^[ (ai + al) sin(2J(,^)) 

+ 2a2(ai - 1) cos(2J(^))) + K3{2a2{ai + 4) sm^{2J{^)) 
+ ((6ai - 3)al + al) ) sin(2 J(0) cos(2J(0) + 2aia2(2ai - 1) cos2(2J(0) ' 
+ U2 (2ciai (a2(ai + 1) sin(2 J(0) + 2aj cos(2 J(0)) ^ 

+ 2^3 cos(2J(0) {a2{al + 2ai + aj) sin(2 J(0) (2(ai - l)al + 2a?) cos(2J(0))) 
+ 2^302^3 ((ai + al) sin(2 J(0) + 2a2(ai - 1) cos(2J(0)' 
- 2w4a2/t3 f 02(^1 + 1) sin(2 J(0) + 2al cos(2 J(0)' 



2«:3a2(ai - 1) (02 sm(2 J(0) + 2a^ cos(2(2 J(0))) ^ 



+ C5, 



+ {{al - ai) sin(2 J(0) + 2aia2 cos(2 J(0)) ^ 

(51) 

where we used Ki = ai + a2, K2 = Oi — 02, ^3 = 0^ + 02- The functions (^j being arbitrary, 
they could be chosen in a way such that their derivatives dQ/d^ simplify the integration 
procedure for F and G. Finally, the solutions for u and v are 

u{x,y) = f{x,y) + F{aix + a2y), v{x,y) = g{x,y) + G{aix + a2y), (52) 

where the functions / and g are given by (H^ while F and G are given respectively by 
the quadratures flSUl) and (15T]) . So, we get a solution of the system ([T]) by defining the 
angle 9 by (fTSj) . (fT9|) . the mean pressure a by (120|) and the velocities u and v by (152|) . 
with ^ (x, y) = aix + 02?/. 

2 Multiplicative separation for the velocities u and w. 

Considering the solution of (IHa), (IHb) for the angle 6 provided by fll5p . 0161) . fll9p . and 
for the pressure cr by (120|) . we now propose the solution for m and to the equations ([He), 
([T]d), in the multiplicative separated form 

uix, y) = f{x, y)F{^{x, y)), v{x, y) = g{x, y)G{^{x, y)), (53) 
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where we suppose / 7^ and (? 7^ 0. For classifying the admissible / and g, we substitute 
the velocities (153|) into the system ([Uc) and ([Hd). This gives 

(a) ifyFi^x, y)) + a,fF\ax, y)) + y)) + a.gG'ia^, y))) sin 2 J(e(x, y)) + 

+ + ai/F'(e(x,i/)) - gyG{0 - a^gG'iO) cos2J(0, 

(&) /.F(e(a;, y)) + ai/F'(e(x, y)) + y)) + a^gCX^^, y)) = 0. 

(54) 

Applying to (l54l b) the linear differential operator defined by (12^ . we have that the 
following equation must be satisfied 



(55) 



where we have F = F{^{x,y)), G = G{C,{x,y)). In dealing with equation (155|) . there are 
three different cases to consider: 

(u) L^{j)=^^ ^Kf)''^' ^^^^ 

(m) ^ 0. 

(i) In the first case, the functions / and g must take the form 

f{x,y)=e'^^^UO, g{x,y) = e'^^^(,{0, ^i^^, (57) 

where Ui corresponds to a damping factor when Ui < 0. Introducing (157|) in the system 
(jSlj), we reduce this one to an ODE system for the functions F and G in term of ^ from 
which we find 

G(e) = Ai(OF(0 + A2(OnO, (58) 
where the coefficients Ai{C,), i = 1,2, are given by 



Alio = 777T -2c^i«2 + cot(2J(0) ~ ( f^T ' 

^2(0 = ^^(Acot(2J(0)-/i), 



(59) 



where (^1, ^2 are arbitrary functions of one variable, while the function F must satisfy the 
ODE 

«^Ci(0^2(0 ^ 

, aMO + a2C( (0^2(0 + a2Ci(0(^i(0 + ^^2(0) . 
+ a2Ci(0^2(0 "^^^^ = 
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If we find the solution for F of the equation (160|) for a particular choice of the functions 
Ci and (2, then the solutions for the system ([T]c), ([T]d), for m and v will be given by fl53l) 
with / and g defined by fl57j) and G by fl58|) . However, by imposing = exp(fr(,^)), 

the problem of solving (160|) is reduced to the one of solving the Ricatti equation 



(61) 



where 

P 8aia2JXO^ _^ 2^(0 ^ 4a;i (ai sin(2J(0) - as cos(2J(0)) J\0 



Cisin(2J(0) C2(0 



Cl 



R rn =^ + (2aia2 J^O + ^1^1 sin(2 J(0)^ - c^^g^ sin(2 J(0) cos(2 J(0)) (^(Q 

'^^^ C2(0 Cisin(2J(0)C2(0 
2a;iJXO(4a2J^(0+u;isin(2J(0)^) 

cisin(2J(0) 

If we choose the function (,2 to be a particular solution of the Ricatti equation obtained by 
imposing that the coefficient -82(0 = then the ODE (IGTD becomes a nonlinear Bernoulli 
equation for K'{^) and the solution for K{^) is obtained by quadrature. In this case, the 
solution for F takes the form 

(62) 



+ £5^2(0- 

(63) 



F(0 = C4 y exp j BiiOd^j + C5, C4, C5 G R. 
We substitute in (155]) and calculate the following solution for G: 



Ai(0 j exp y 5i(0t^e) rfe + ^(Oexp y Si(e)c^e 



G(0 = C4 



Finally, the velocities u and f are 



C2{y/x) C4 / exp 



u(a;,?/) 



v{x,y) = e--'^Uyl^)\cAAr{y/x) 



^=y/x 



5i(0rfnrfe + C5 



(64) 

where Ai, ^2 are defined by (159|) . Ci(0 arbitrary function and C,2{0 ^ solution of 
the Ricatti equation 

4J'(0 (2aia2 J'(0 + wiOi sin(2J(0)' - ^102 sin(2 J(0) cos(2J(0)) ^(0 



C^'(0 + 



Cl sin(2J(0) 



^ 2c^iJX0 (4a2J-(0+a;isin(2J(0)^) ^ 
Cl sm(2J(0) 

So, we have obtained a solution of the system ([1]) by redefining the angle 9 by (fT5|) , ([in 
the pressure cr by (120|) and the velocities u and f by with ^^^^ 1/) = c^i^^ + 
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(ii) In the second case, we have that 

g{x,y) = uj,e^-'^Ci{0. c^3,c^4eR, (65) 
and the function f{x,y) must be a solution of the system of two PDE 

/. = -aiC2(0/ + ^7C3(0 

-2a2 cos(2J(0) + ai sin(2 J(0) 



(66) 



aa sin(2J(0) 



where the Ui, i = 1,2,3,4, are real parameters and, in order that the system (166|) be 
compatible, the functions 0(0; ^ = 1)2,3, must obey the condition 

2Ci(OC3(0 { ^2i2j[f)) - ^"^^ sin(2J(0) + a2Cos(2J(0)) 
+ (Asin(2J(0) + /icos(2J(0)) (Ci(0C3(0 + diOUO^O + Ci(OC((0) (67) 

uJi \ a2 J J 

The functions F and G are given by 

F(0 = u;ie^^«), G(0 = ^ ^ Ci(0C3(e)e^^«)rfe) • (68) 

As an illustration, we compute an explicit solution of the system (fT]). To do this, we 
choose 

C2(0 = -ln(Ci(0), C3(0 = 2a;5cr'sin(2J(0)J'(0, ^5 G R, (69) 

so that the equation flB7|) is verified and consequently the system can be solved 
because it becomes compatible for /. The solution of the system (166|) is then 



(70) 



where ^ = aix+a2l/ and the function Ci(0 has as its only restriction that its first derivative 

Finally, by introducing (153]) . fl68p and (170]) in the expression flS^ . with the function 
6'(a;, y) defined by f[T^ . (ITBl) and f[T^ . the solutions for m and v of the system ([He), ([Hd), 
are 

M(a;, y) = -C4a2 cos(2 J (0) + CiC4y + Cg, 

f (x, y) = c^ai cos(2J(^)) + ce, q G R, "i = 1, 2, 3, 4, 
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where the real parameters Ui have been redefined to simphfy the expression flTTl) in which 
the new parameters are q. So, we have a solution of the system ([1]) by defining the 
angle 6 by (fT5|) . (fT9|) . the pressure cr by (120|) and the velocities m and f by fITT]) . with 

i{x,y) = aix + a2y. 

(iii) If the condition (13^ iii) is satisfied, then we have 

F{0 = uj,e-^'^^y'^\ G(0 = -c^lC2(0e"^^^^^/''^ e R, (72) 
when the functions f{x,y) and g{x,y) are solutions of the system 

fy =-C((0/ + C2(0^7. - 2cot(2J(0)C2(0^7y (73) 
ai 

- (aia2-i-2cot(2J(0))C3(ek, 
where the functions Ci(05 C2(0 arbitrary and 

C3(0 = a2 (- "-'^^If^K my (74) 

Using the compatibility condition on the mixed derivatives of the function f{x, y) relative 
to X and y, we obtain the second order ODE for g{x,y) 

C2(0 {9xx - 9yy - 2 COt (2 J(0 ) <?.,) - 2(3(0 («1«2"' " COt(2J(0)) 

+ 2 (-C3(0 (l + aia,-^cot(2J(0)) + ^^'^'^^^ 



sin^(2J(0) 



+ 



sm (2J(4)) aia2 J \ 0.2 J 



(75) 

It's a hyperbolic equation on the domain Vl defined by fl^T]) . So the change of variable 

_ 2a2J(0 , ai(cos(2J(0) + l) + a2sin(2J(0) 
[aj (p[x, y ) — X ^ -\- ^ , 

, 2aJ(0 ^ ai(cos(2J(0) + l) + a,sin(2J(0) ^ ^^^^ 

Ci Ci 

transforms the equation (175|) to the simplified form 



Ci / COS(2J(0)-1 , C3(0 ^ 

= ^.,(0, ^) + - ,i,(2J(0) + a.C2(0^'(0 J 

^ / cos(2J(0) + 1 _ ^3(0_\ . , 
+ 4aA sin(2J(0) a.^^AO J ^''^^'^^ 

ci aiC3(0 , Asin(2J(0) -/icos(2J(0) i i\ 

+ A„ „ ^ Tuc\ (C3(0Ci(0 - 01(3(0) ) 9{<f), ^) 

(77) 



2al VC2(0sin(2J(0) 4a2aiC2(0^'(0 
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Once again, if we take the difference of equations fl76la) and fl76lb). we deduce tfie relation 
(H5|) . In order to simplify the equation (177|) and to solve it, we choose C3(0 — 0- This is 
equivalent, considering (17il) . to requiring — a]~^C2(0Ci(0 + (2(0 = O5 which is satisfied if 
we choose 

C2(0 = ^3exp(ar^Ci(0)- (78) 
So the equation (177|) transforms to 

9<P^ + , "'w't,!^ - f «i sin(2 J(</>, ^)) + :^(a^ - a?) cos(2 J(</., ^)) 
4a2sm(2J(0, ■i/')) \ 2a2 

ca A ^^^^ 

^ sin(2 J) (-A sin(2 J) + /i cos(2 J(0, V^))) j " ^'^^ " °' 

We propose the solution of (179|) in the form (7 = g{J{(f), ip), K{(j), tp)) with J given by (H5|) 
and = (2a2)~^Ci(^/' + 0), which allows us to find the solution 

g{<P, ^) = Us cos(2 J(0, ^)) + sin(2 J(0, ^)) - ^^^^^ + V^) ^ ^g^^ 

The solution expressed in term of x and y is 

. X / ^^4^1 ^ /o7^c^^ Uijxci + ai) 
g{x,y)=[uj3 cos(2J(0) + ws- (81) 

\ 0-2 / 0-2 

Since the compatibility condition fl7^ is satisfied, the system (175|) can be integrated after 
the substitution of g{x,y) given by (1ST]) . C2(0 by (175]) and C3(0 = 0. The solution for 
f{x,y) is given by 



= gCi(0 ( ^3(-Q2C^3 + ^4ai)cos(2J(0) _ Cia;gy ^ 



2„ 

82) 

ai ai y ~ 

where Ci(0 arbitrary function of one variable and C, = cliX + a2y. Finally, the 

solution for the system ([He), ([Hd), is obtained by introducing (172]) . (IHT]) and (182]) in (153]) . 
Redefining properly the parameters Ui, i = 3, 4, 5, 6, the solutions are 

u{x, y) = C4 (-as cos(2 J (0) + c^y) + Cg, ^^^^ 
f (a;, y) = C4 (ai cos(2 J(^)) + (-C5 + ci)x) + C7, 

where the Ci, i = 1, . . . ,7 are integration constants. So, we have obtained a solution of the 
system ([1]) by defining the angle 6 by (|T5]) . (|T9]) . the pressure a by (120]) and the velocities u 
and f by (15^ . with ^{x, y) = aix + a2y. This solution includes the solution (I7ip to which 
it corresponds if we choose C5 = Ci in (155]) . A extrusion die corresponding to the solution 
( 155]) is shown in figure El The parameter values that have been used are k = 0.8, Ci = 2, 
C3 = 1, Oi = —10, 02 = 10, C4 = 1, C5 = 1, Cg = 10, C7 = —10. To demonstrate how the 
plasticity region changes when we vary the feeding velocities, we plotted two curves Ci 
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and C2 to delimit the plasticity region at the mouth of the tool. The curve Ci corresponds 
to a feeding velocity Uo = 3.54 and the curve C2 to a feeding velocity Uo = 3.84, while the 
curve C3 is the plasticity limit at the output and the extraction velocity is Ui = 16.12. 
This kind of tool could be used to thin a plate for some ideal plastic material. 



Figure 3: Extrusion die corresponding to the solution (PO)) . (|83p . 



B Similarity solution for the angle 6 and corresponding pressure a 

In this section, we find solutions of the system ([T]) for which the angle 6* is a similarity 
solution. We propose the solution for 6 in the form provided by (llSp . but with the 
symmetry variable in the form 

^{x,y) = y/x. (84) 

The solutions for 6 and a, which are obtained by this Ansatz, include for appropriate 
parameters choices the invariant solution corresponding to the subalgebras Cij, j ranging 
from 5 to 10 and L2,k for k ranging from 13 to 27. The introduction of f llSp . with ^ defined 
by (I84p ■ in the system (fHa), ([Hb), leads to the system 

a^{x,y) = 2{k/x) {-^{x,y) cos{J{^{x,y))) + sm{^{x,y))) J'{^{x,y)), ^^^^ 
(^y{x,y) = -2{k/x) {^{x,y) sm{J{^{x,y))) + cos{J{^{x,y))))J'{^{x,y)). 

Considering the compatibility condition on mixed derivatives of a relative to x and y, we 
deduce from fl85|) the following ODE for the function J: 

ae-l) sin(2 J(0) + 2e cos(2 J(0)) J"{0 + 2 sin(2 J(0) + cos( J(0)) J'iO 
+ 2(-2sin(2J(0) + (e'-l)cos(2J(0))(J'(0)' = 0, ^ ^ 

which has the first integral 

ae - 1) sin(2 J(0) + 2^ cos(2 J(0)) J'(0 = c^, (87) 
where Ci is an integration constant. There are two cases to consider to solve the equation 

dHZD. 



i If Ci 7^ the solution of (187|) is given in implicit form by 



(tan(J(O)-Ov/^ 
(tan(J(0)e + l)ci-e + tan(J(0) 
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where C2 is an integration constant. The solution for a is obtained by integrating the 
system (185|) and taking into account the first integral (187|) . We find 

y) = k y) cos(2 J(e(x, y))) - sin(2 J(e(x, y))) - 2ci ln(a;)) + C3, C3 G R (89) 

with ^{x,y) defined by (15^ and J(^) a solution of ( ISSj) . So, 9{x,y) = J{y/x) and <j{x,y) 
are given by (I89p and are solutions of the system ([Ha), (IHb), if the function J(^) verify 
the algebraic equation 



J(^) = - - arctan ( ) . (90) 



ii If the constant Ci = in the equation fl88p . then the solution of (1861) for J is 

1 arctan (^ 

We subsequently solve (1851) for a{x, y) considering (I90p . which leads to the solution of the 
system ([Ua), ([Hb), given by 

1 f 2,xy \ fy\ 

6{x, y) = - arctan I — I , a{x, y) = —2k arctan ( — ) + C2. (91) 

2 \x^ — y^ J \x/ 

1 Additive separation for the velocities when ci ^ 0. 

Already knowing the solution 6{x,y), a{x,y), in the case of ci 7^ 0, we still have to 
compute the solution for u and v. A way to proceed is to suppose that the solution is in 
the additive separated form 

u{x,y) = f{x,y) + F{^{x,y)), v{x,y) = g{x,y) + G{^{x,y)), (92) 

where C,{x,y) = y/x. We introduce (192|) into the system ([He), ([Ud), which gives 

(sin(2 J(0) - cos(2 J(0)) ^F'iO - (cos(2 J(0) + sin(2 J(0)) ^G\0 .gg. 
+ ((/, + ^7,)sin(2J(0) + (/.-^7,)cos(2J(0))x = 0, ^ ^ 

(/x + ^7,)x + G"(O-ei^'(O = 0, (94) 
We must determine which functions / and g will reduce the equations (1930 . (I94p . to a 
system of ODE for the one variable functions F{C,) and G{^). To reach this goal, we 
first use as annihilator the infinitesimal generator Di defined by ([3]) that we apply to 
the equations (193|) . flMj) . to eliminate the presence of the functions F and G. Indeed, 
Di annihilates any function of ^ = y/x. So, we obtain as differential consequences some 
conditions on the functions / and g. We can assume fx{x,y) + gy{x,y) 7^ 0, otherwise 
we can show that the only possible solution is the trivial constant solution for u and v. 
Under this hypothesis, the previous conditions read 

f. = -gy + Ci{Ox-\ (95) 

fy = -9. + 9yC2{0 + (96) 

{gy + D, (gy)) (C2(0 sin(2 J(0) - 2 cos(2 J(0)) = 0, (97) 
where the functions of one variable Q, i = 1,2,3 are arbitrary. The left member of ([9j 
being composed of two factors, we must consider two possibilities. 
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a. We first suppose that 

Qy + {gy) = 0. (98) 

In this case, we find that the functions / and g take the form 

f{x,y)=-J + u^\n{y} - uiy + U5, ^gg^ 

g{x, y) =(4(0 + ^2 ln(x) + uix + W3, 

where the functions C,i{C)iCi{C) arbitrary and the functions C2(0? CslO? were chosen 
to solve the compatibihty conditions on mixed derivatives of / relative to x and y. We 
now introduce the solution (199|) in the system (l95l) . (196|) . which leads to an ODE system 
for F and G, that we omit due to its complexity, and for which the solutions take the 
form of quadratures 



FiC)-f( C'(n I ^'^^^ I (c^4 + c..0sin(2J(0) \ 

y V ^'^^^ e +e((e-l)sin(2J(0) + 2ecos(2J(0))r^ + '^' 



^(0 



^ y (-Ci(0 + ei^'(O) + C4, C4, C5 e R. 



(100) 



The last step is to introduce (199|) and fllOOp in the Ansatz fl92|) . Then the velocities are 

l c5Cos(2J(y/x)) , /•^/- C6sin(2J(0)JXO ^^ , 
M(a;, y) = - + / + cg \n[y) - ay + C7, 

v{x,y) =C5ln(x) +C4X + ^^^cs^ sin(2 J(0) /(O^^ - ^^^^^^^^ + cg/, 

where the q are integration constants. So, we obtain a solution of the system ([1]) by 
defining the angle 9 by ffTSj) . fl88|) . the mean pressure o" by (189|) and the velocities u and 
f by fllOip . with ^(x, y) = y/x. 



b. Suppose now that the condition (l97j) is satisfied by requiring 

C2(0 = 2cot(2J(0). (102) 

Then, applying the compatibility condition on mixed derivative of / relative to x and y to 
the equations (195|) . (196|) . and considering (^2 given by (I102p . we conclude that the function 
g must solve the equation 

gU^, y) + 2 cot(2 J(0)(7.,;(x, y) - gyy{x, y) + 4^"^^ [i + cot2(2 J(0)) /(0^7y(a;, 1/) 

+x-\c',m+m)+m) = ^- 

(103) 
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It's a hyperbolic equation everywhere in the domain where J is defined. So, we introduce 
the change of variable 

. N , /'^^"'''^ sin(2J(0) 

1 + cos(2 J(4)) + 4 sin(2 J(4)) 

(104) 

^ , f^^"''^ sin(2J(0) 

-1 + cos(2J (4)) + 4 sin(2J (4)) 

which brings the equation fll03p to the simplified form 

ci / sin(2J(0,^))5(0 sm{2 J {(t>,tp))g^ 



''^ ' 2 VV'(cos(2J(0,^/'))) + 1 0(cos(2J(0,^/')) - i; 
1 sin(2 J(0, ij))e - sin(2 J(0, ^)) + 2^ cos(2 J(0, ^)) ((3(0 + C[ (0 + iCm) 



0, 



4 

(105) 

where J{(j),tp) is defined by fl45p . To solve the equation fllOSp more easily, we define the 
function (3 by 

C3K) . + (106) 

So, the solution of fllOSp is 

^(0, ^) = -(1/2) (a;i - (1/2) ln(^/</))) cos (ci ln(7/>/7/>)) - (l/4)cr' sin (ci ln(^/0)) , (107) 
which, returning to the initial variables, takes the form 

g{x,y) = -(1/2) {ioi - (l/2)q V(i//x)) cos(2 J(y/x))-(l/4)crSin(2J(t//x))+u;2. (108) 
After the introduction of the solution (11081) for g, the function / is given by quadrature 



from the equations (l95j) . (196|) . The obtained solution for / is 

f( , , /•^^^•^ Ucia;i-J(0)sin(2J(0)JX0 ,A /"^^^''^ , w 

/(x,?/)=| / — ^)~J —^d^+{ci+l)uJi\n{y)+uj3. 

(109) 

We now introduce fll08p . f llOQp . in f l93|) . (IMj) . and get F and G by quadrature in the form 



^ " J ci^ " (110) 

G(0 = -(l/2)q ^ ((u;i + J(0) cos(2 J(e)) + 2 sin(2 J(0)) • 

Finally, the substitution of ffTOSD . ffT09|) and ffTTOj) in ([92]) provides the solution to ([Uc), 
dUd): 

. N . ,^ w N r/"c4(ci + l)sin(2J(0)J'(0 
m(x, y) =(ci + 1)C4 ln(|/) + y L_}_^^llL^d^ + C5, 



C4(ci + l)cos(2J(?//x)) 
^(a;,!/) = h C6, 



111 
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where the q are integration constants. So, we have a solution of the system ([T]) by 
imphcitly defining the angle 9 by (fT5|) . flHHj) . the mean pressure by a by (189|) and the 
velocities m and v by fillip , with ^(a;, = y/x. 



2 Additive separation for the velocities u and v when ci — 0. 



Now, we consider the case where Ci = in (l87j) . Then the solutions for 6* and a are given 
by (19T|) . We still suppose the solution for m and v in the form (192|) . The procedure is the 
same as the previous case until we obtain the conditions (195|) . fl96|) and (197|) . We must 
again consider two distinct cases. 



(a.) We first suppose that the condition (l98l) is satisfied. Then, the functions / and g 
are defined by 

gix, y) =(4(^(3;, y)) + ujix + UJ2, e R, i = l, 2, 3, 

where the (^j, = 1, 2, 3, are arbitrary functions of one variable and to simplify the expression 
for / and g we have chosen Ci(0 = (1 - (2(0^ - C4(0 -^CalO- We substitute (fTT2D in 
the equations (193|) and to determine F and G. We conclude that F{^) is an arbitrary 
function, while G{S,) is expressed in term of a quadrature by 

G{i) = J {F'm + (ec2(o -i+e)c- ec3(o) (ns) 

We finally obtain the solution u and v by introducing (I112p , (I113p in ( l92|) and, to simplify, 
by choosing 

C3(0 = -(C2(0 + OC(0, 

which gives 

u{x,y) =-C2y + F'{y/x), 

v{x, y) =C2X + {y/x)F'{y/x) - F{y/x), 

where F is a arbitrary function of one variable. A solution of the system (fl]) consists of 
the angle 6 and the pressure a defined by (1911) with the velocities defined by (11141) . For 
example, if we choose the arbitrary function to be an elliptic function, that is 



F(0 = dn (^biTT (1 - , g 



and we set the parameters as k = 0.027, e = 1, = 0, 61 = 2, 62 = 2, ^ = 1/2, then we 
can trace (see figure Hj) an extrusion die for a feeding speed of Uq = 0.75 and an extraction 
speed Ui = 1. The curve Ci on the figure H] delimits the plasticity region at the mouth 
of the tool, while the a;-axis does the same for the output of the tool. This type of tool 
could be used to undulate a plate. We can shape the tool by varying the parameters. 
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For example, we can vary the wave frequency setting the parameter bi. Moreover, one 
should note that if the modulus g of the elliptic function is such that < < 1, then 
the solution has one purely real and one purely imaginary period. For a real argument 
we have the relations 

Vl - p2 < dn{x, q) < 1. 




Figure 4: Extrusion die corresponding to the solution (QUI) . (|114p . 



b. Suppose that C2(0 defined by fll02p and for simplification we choose in particular 

C3(0 = (115) 



Applying the mixed derivatives compatibility condition of / to the equations (195|) . (196|) . 
we get the following ODE for the function g: 

By the change of variable 

i{.x,y) = y/x, T]{x,y) = x'^ + y'^, (117) 

we reduce the PDE flll6p in term of x and y, to the much simpler PDE in term of ^ and 
V, 

9,, + ^ = 0, (118) 

which has the solution 

9iC.v) = U0 + 4==, (119) 

+ 1 

where (4 and are arbitrary functions of one variable. Then, we find the solution for / 
by integration of the PDE dHS]), ([96]), with Cs given by ffTT5|) . 

f{x,y) = - T y , , +C2. (120) 

J ? vv{x,y) 
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By the substitution of flllQp and fll20p in the equations 
arbitrary function of one variable and G is defined by 



931), (El, we find that F is an 



:i21) 



So, we introduce (IllOp . (I120p and (I12ip in ( 192|) and after an appropriate redefining of (i, 
(4 and (5, the solution to ([He), ([Hd) is provided by 

u{x,y) = K'{y/x)-yH{x^ + y^) + C2, v{x,y) = -K\y/x)+^K{^)+xH{x^ + y^) + ci, 

(122) 

where H, K are arbitrary functions of one variable. The velocities (I122p together with the 
angle and pressure defined by (19T]) solve the system ([1]). For example, a tool corresponding 
to the solution ( I122p with H{ri) = 2 exp(— O.lr^), K{C,) = ^ and for feeding and extraction 
speed given respectively by Uq = 1.05 and Ui = 1.05. It is shown in figure O The 
plasticity region limits correspond to the curves Ci and C2. This tool is symmetric under 
the reflection x 1— )■ —x. Moreover, the top wall of the tool almost makes a complete loop, 
and this lets one suppose that we could make a ring in a material by extrusion. 




Figure 5: Extrusion die corresponding to the solution ^U\i. p22p . 



3 Multiplicative separation for the velocities u and v when ci 7^ 0. 

Consider the solutions of ((Ha), ([Hb), given by the angle 6{x,y) = J{y/x) with J defined 
by (188|) and the pressure a defined by (189|) . We require the solutions for u and v in the 
multiplicative separated form 

uix, y) = f{x, y)F{^{x, y)), v{x, y) = g{x, y)G{^{x, y)), (123) 

where ^{x,y) = y/x and /, g, F, G are to be determined. The Ansatz (I123p on the 
velocities brings the system ([Uc), ([Hd) , to the form 

[fyF + x-'fF'iO + 9MO - x-'^gG'iO] sin(2 J(0) 

+ [f^FiO - x-^UF' - gyG - x-^gG'] cos(2J(0) = 0, (124) 

[/.F(0 + gyG{i)] - UF\i) + gG\i) = 0. 
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To reduce the PDE system fll24p to an ODE system involving F and G in term of ^, we 
act with the operator Di, defined by ([3]), and annihilate the function of ^ present in fll24p . 
This leads to conditions on / and g that do not involve F, G and their derivatives. There 
are three cases to consider, that is 

(a) D^{gy/g)^Q, 

ih) Diixgy/g)^0, 
(c) Di{xgy/g) = 0. 
In this paper, we present the details for the cases (a) and (c). 



Di 



Diif/g) 

Di{xgy/g) 

Diif/g) 
Di{xgy/g) 



7^0, 



(125) 



(a) We first suppose that the condition fll25[ a) is satisfied. In this case, the function / 
must be solution of the PDE system 

/x = - x-'^CiiOf + C2i09y - (C2(0CU0 - cm) 9, 
fy =x''Ci{Of + C2{09. - 2C2(0 cot(2J(0)^7, .^26) 
-u (r- m^'m ^'m^ esin(2J(0) + 2cos(2J(0) 

where (i, (2 are two arbitrary functions of one variable. For the system fll26p to be 
compatible, the function g must satisfy the PDE 

g,, - 2 cot(2 J(0)^?., - gyy + 2x-iC2(0"' {i + cot(2 J(0)) UO9. 
+ 2x'^ [{i - cot(2J(0)) C3(0 - 2esin-2(J(0) J'(0] 9y 

+ x-\2{i)-' {miO sm-^(2 J(0) + - 1 + 2e cot(2 J(e))) C( (0 ^^^^^ 
-2^-2 cot(2 J(0)]C3(0 - - 1 + 2^ cot(2 J(0)) CslO)^? = 0, 



where we used the notation CslO — C2(0Ci(0 ~ C'li'O shorten the expression. The 
equation (I127p is difficult to solve for arbitrary (^1, (,2-, but if we make the particular choice 

C2(0 = ^^2exp(Ci(0), c^2eR, (128) 

then the PDE ffT27D reduces to 

g,, - 2 cot(2 J(O)^?.,, - gyy - ^x-^iJ'{i) sm-\j{i))gy = 0, (129) 

which is solved by the function 

5'(x, y) = + (X'4 cos(2 J(?//x)) + Ws, Wa, u;4, ^5 G R. (130) 
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With g given by fll30p and C2(0 by fll28p . the system fll26p is compatible. Consequently, 
/ is expressed in term of a quadrature. We find 

f{x,y)=2exp{Uy/x))uj2Uj, j'' M'^m)J'iOd^ (^31) 

+ exp(Ci(?//x)) (^2^3 + 2uj2UJiCi ln{y) + Uq) 

With / given by fllSip and g by fll30p . the solutions for F and G of the system fll24p are 

F(0 = cui exp(-Ci(0), GiO = C02. (132) 

By introducing fll30p . fll3ip . f ll32p in f ll23p and redefining the free parameters Ui, i = 
1,2,3,4, the solutions of ([He), ([Hd), for the velocities u and v is 

/y/x 
r' sin(2 J(0) J'(Oc?e + C5y + 2ciC4 \n{y) + cg, ^^^^^ 

v{x, y) = —C5X — cos(2J(y/x)) + C7, 

where the q, i = 1, . . . , 7, are integration constants and J is defined by (IHHl) . So, we 
have a solution of the system ([1]) composed of the angle 6 in the form (fT5|) with J given 
implicitly by (188|) together with the pressure a ( 189|) and the velocities f ll33p . 

(c) Suppose now that the condition ([1251 c) is satisfied. In this case, the solution g takes 
the form 

g{x,y) = h,{x)CM^,y)) (134) 
and the function / must be solution of 

fx=X-'^C2iOf + ^'\3mii^), 

fy = - 2x-i/ii(x)C3(0 cot(2J(0) + e-^^(«) ^ e^^«)C3(0^e - h[{x) (135) 

-x-^ [C^(0 + eC3(0]- 

We omit, due to its complexity, the expression of the compatibility condition on the mixed 
derivative of / relative to x and y. Nevertheless, making the specific choice 

the system (I135p turns out to be compatible and the solution to (I135p is 

fix, y) =Uylx) f (1 - [y/xf - 2{ylx) cot(2 J(y/a;))) f'^ sin(2 J(0)/(0^^^ 



(137) 



+ [1 - {y/x^ - 2{y/x) cot{2 J{y/x))] {\ny - 00^] 
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The substitution of fll37p . f ll34p and fll36p in the system fll24p results in an ODE system 
for F and G, omitted due to its complexity, which has the solution 

p^^^ u;isin(2J(0) 



C3(0 ((e - 1) sin(2J(0) + 2^cos(J(0)) ' (138) 
G'(0 = Ci(0'M^icos(2J(0)-u;2). 

We finally, obtain a solution for the system ([He), ([Ud), by the substitution of fll37p and 
ffT38|) . with hi{x), (2(0. defined by ([^6]) . in ffT23p . This leads to 

u{x, y) = -C4 \n{y) - c{\^ J'^ sm{2J{0)J'{0d^ + C5, ^^gg^ 

v{x,y) = (l/2)q ^u;icos(2J(y/x)) + Cg, ci, 04,05,05 G R. 

So, the system ([1]) is solved by the angle 6 in the form (fT5|) with J implicitly defined by 
together with the pressure a fl89|) and the velocities fll39p . 



4 Multiplicative separation for the velocities u and v when ci = 0. 

Consider now the case where Ci = in (157|) so the solution of ([Ha), ([Hb), for 6 and cr is 
(I9T]) . We suppose that the velocities m and f are in the form fll23p . Introducing this form 
for the velocities and 9 defined by (19T]) in the equations ([He), ([IJd), leads to the system 
fll24p which reduces to a ODE system for F and G if the functions / and g satisfy the 
condition fll25p . The three different constraints fll25p must be considered separately. 



(a) In the first case, where we consider that the conditions fll25[ a) are satisfied, the 
functions / and g must verify the system (11261) . (I127p . Changing the variables {x,y) to 
new variables {C,, if) defined by 

^ix,y) = y/x, r]{x,y) = x"^ + y'^, (140) 

and considering 6 given by (1911) . the system (I126p . (ll27p . becomes 

= -Ci(o/ - r'C2^5 - c ((2(0 - Ci(0C2(0) 9, 
— + + 1)-' = 0, (142) 

9rj (2(0 
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where Ci(0? C2(0 arbitrary functions of one variable. The solution of the system 
fll4ip . f ll42p . for / and g as function of ^ and rj is 

f{i,v) HUO {K{i) + Ui)-\^ + e)-"' (//(r/)e-/^^«)'^«)) 

_e'/^^«)'^q j r^(e2 + i)[c2(o + (Ci(o+e(e' + 1)^^(2(0)^(0 

^ (143) 

+ C2(0^'(0]rfe-^ij, 

^7(e,r/) =K{i) + (e/^^«)'^«C2(0v/e+T) , 

where the functions K{^) and -/^(r/) are arbitrary functions of one variable. Introducing 
/ and g expressed in the initial variables y, by the substitution of (I140|) in (I134p . the 
system (I124p is reduced to an ODE system for the functions F and G in term of ^ which 
have the solution 

F(0 = ooie^^'^^^''^ G(0 = -cJiC2(0e^^'^^^''^- (144) 

Finally, redefining 



KiOr' exp (^J QiO, H{r]) = co^'^Piv) 



C2(0 = -(^1 

where Q{^), P{ti) are arbitrary functions, and doing the substitution of fll43p and fll44p 
in fll23p . we obtain the solution of ([He), ([Hd), for velocities u and v given by 

uix, y) = yPix' + y') - x-'yQiy/x) j''^ ((^^ + 1)Q'(0 + ^^(0) + ^2, ^^^^^ 
v{x,y) = Q{y/x) - xP{x^ + y'^), 
This solution for velocities together with 9 and a defined by (191]) solves the initial system 

m- 

(b) Consider now that the functions / and g satisfy the constraint (11251 b). Then they 
take the form 

f{x,y) = Ci{ylx), g{x,y) = C2{y/x), (146) 

where (i and ^2 are arbitrary functions of one variable. Since the functions / and g 
depend only on the symmetry variable and the velocities u and v have the form (I123p . it 
is equivalent to consider that 

u{x, y) = F(Ox, y)), v{x, y) = Gi^^, y)). (147) 
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If we suppose that u and v are in the form fll47p . then the solution of ([T]) consists of the 
angle 6 and the pressure a given by (19T|) together with the velocities 

uix, y) = Fiy/x), y) = j'' iF\i)di + cg, (148) 

where F is an arbitrary function of one variable. 



(c) The third case to consider is when / and g obey the conditions fll25[ c). so they take 
the form 

f{x,y) =uj,x'^--{y/x)^'+--'^'\x-\x^ + y^)r'\ g{x,y) = x'^^^Uv /x). (149) 

where Ci(0 arbitrary function of one variable. Then we introduce (I149p in (I124p 

and solve for F and G. The solution is 

F(o = ujic^\m^'-^'^'\e + ^(o = {e + ir"- {m 

Finally, substitution of /, F, G, given by fll49p and fllSOp . in fll24p gives, after redefining 
the parameters Ui in a convenient way, the solution for u and v of the equations ([He), 
dHd), 

uix, y) = c^y{x^ + y^^'^ v{x, y) = -c^x[x^ + y'Y^I'' , (151) 

where C3, C4 are integration constants. The velocities u and v together with the angle Q 
and the pressure a given by (19T|) constitute a solution for the system ([1]). This solution is 
just a subcase of the previous one corresponding to the condition (1125 l a) and the choice 
Q(^) = 0, P(r/) = r/'^2/2 in ([T15D . 



C Solution for in terms of the invariant a{x,y) + aix + 02?/. 

The subalgebras £2,28 to £2,31 do not have symmetry variables, that is an invariant de- 
pending on X and y only, but they all have an invariant in the form a + aix + a2y for 
appropriate values of ai, 02. So, we suggest 6 in the form 

eix,y) = JiT{x,y)), (152) 

where 

T{x,y) = a{x,y) + aix + a2y. (153) 

The function t{x, y) is to be determined since it depends on one of the unknown quantities, 
that is a{x,y). Consequently, finding the function r also determines the pressure a{x,y). 
We first introduce 6 given by (I152p in the system ([Ua), ([Ub), which leads, after the 
elimination of a{x,y) with the use of (I153p . to the system 

ai + 2A; (ai cos(2J(r(x, y))) + 02 sin(2J(r(x, y)))) J'(r(a;, I/)) 
^ l-4PJ'(r(x,y)) ' 

. . _ a2 + 2A; (ai sin(2 J(r(x, y))) - 02 cos(2 J(r(x, y)))) J'{t{x, y)) 
""^^"^'^^ ~ l-4A;2J'(r(x,y)) ' 
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of the two PDE for the function t{x, y) using the compatibihty conditions on mixed 
derivatives of r relative to x and y, we deduce the second order ODE for the function 
J(r) 

4k {al + al) 



^ ^ (ai - «?) sin(2 J(r)) - 2aia2 cos(2 J(r)) ^""^^ 
2aia2sin(2 J(r)) + {aj - a^) cos(2 J(r)) ax2 ^ ^ 
(ai - af ) sin(2 J(r)) - 2aia2 cos(2 J(r)) ^ ^^^^ 
The solution of the previous equation is 



(155) 



J(r) = i arctan ( ^^^^3^^ ^ ) , (156) 

where 

A = (ci/2) — af) , fx = k{a\ + a\){a\ — a\)^^ — Ciaia2, ci G R. (157) 
The solution of (I154p for r takes the implicit form 

/ti sin(2 J(r + C2)) + fi;2 cos(2 J(r + C2)) + X + + C3 = 0, 02,03 eR, (158) 
where 

2kai c^Aa^ — alY — Ak'^ 0102(0? — a^) + 2A;^ai , , 



2ciai(af — a^) + 4A;a2 ' Ciai(af — aD + 2A;a2 



Then we introduce (I156p in 0158p and we solve for r to obtain the solution in the explicit 
form 



r (X, y) = -^^^±^{x + K,y + C3) ± V^ i + " + +^3)^) (A^i - ^^2)^ ^^^^^ 

where ki,K2, K3 are given by (I159p and A,/i are defined by (I157p . Finally, the solution to 
the PDE (Pa) and (fUb) is 

6'(x,y) = J(r(x,?/)), a(x,?/) = r(x,y) - (aix + 02?/), (161) 

with T{x,y) defined by (I160p and J by (I156p . 

1 Additive separation for the velocities u and v. 

Finally, we search for a solution to the PDE (JHc), (fUd), in the additive separated form 

uix,y) = f{x,y) + F{T{x,y)), v{x,y) = g{x,y) + G{T{x,y)), (162) 

where T{x,y) is given by fll60p . We substitute fll62p in the system ([He), ([Ud), to obtain 

(/, + TyF'ir) +g, + T.G'ir)) sin(2 J(r)) + (/, + t,F'{t) ~ - r,G"(r)) cos(2J(r)) = 0, 
f. + r,F'{T)+gy + ryG\T) = 0. 

(163) 
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We want this system to reduce to a ODE system for r. This is the case if the functions 
/ and g satisfy some conditions. In order to find these conditions, we act on the system 
f ll63p with the annihilator of r defined by 

Lt = -Tyd^x + r^dy. (164) 
The resuh is that / and g must obey the PDE system 

fx = Ur)-gy, 

fy = -g, + 2coi{2J{T))gy + C2{r), ^ ' 

where ^i, C2 are arbitrary functions of one variable. The compatibihty condition on mixed 
derivative of / relative to x and y leads to the following PDE for g: 

g^^ - 2 cot{2 J {T))g„j - gyy + —f^2J^9y + Ci{.r)Ty - C2(r)r^ = 0. (166) 

Any solutions / and g of the system consisting of f ll65p and f ll66p reduce the system fll63p 
to an ODE system for F and G in term of r. The general solution of f ll66p is hard to find, 
but we give as an illustration a particular solution for g. We make the specific choice 

> / N ^ , s 2u;2 cos(2 J(r)) , 

0M=.„ C.W = -— pyy^+-3, (167) 

where Wi, u:2 and are free real parameters. Considering C,i{t), C2(t) given by fll67p . the 
solution for g{x,y) of the equation f ll66p is 

g{x,y) = u^x + ijJ2y, E H. (168) 

We introduce fll67p and f ll68p in the system fll65p and we solve it to find 

f{x,y) = (ui - uj2)x + {-uj^ + uj3)y + U5, UqER. (169) 

Using f ll68p and f ll69p . the solutions F{t) and G{t) of the system fll63p are 

F{t) = K4, sin(2 J(r)) + cos(2 J(r)) + C4, G(r) = kq sin(2 J(r)) + K7 cos(2 J(r)) + C5, 

(170) 

where C4, C5 are integration constants and 

/ 2kai \ a2Ui - aiUs 

= [ui - U2) I a2Ci + , K5 = {l/2){aiUJi - a2UJ3)Ci H ^ 2 — ' 

\ a]^a2 / di CL2 

2ka2 02^3 - a-i^i 
Kg = (aiCi 2_ 2 )^2, = (l/2)(aia;3 - a2Wi)ci H 2 _ 2 — • 

^1 ^2 ^1 ^2 

Finally, the solutions of ([He), ([Hd), for the velocities u and f are 

u{x, y) = {ui - U2)x + (-W4 + U3)y + K4, sin(2 J(r(x, y))) + K5 cos(2 J(r(x, y))) + C4, 
f (x, y) = Co'4X + U2y + Kg sin(2 J(r(x, y))) + K7 cos(2 J(r(x, y)) + C5, 

(171) 

where the function J is defined by fll56p . The velocities fll7ip together with the angle 6 
and the pressure a defined by (I16ip in which J and r are respectively defined by (I156p 
and f ll60p solve the system ([1]). 
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IV Concluding remarks and future outlook 

The main goal of this paper was to construct analytic solutions of the system ([1]) modelling 
the planar flow of a ideal plastic material and to use the flow described by these solutions 
to deduce possible shapes of extrusion dies. The group theoretic language appears to be a 
very useful tool to obtain these solutions in the analysis of their admissible flow and their 
properties. The main benefit of using group analysis is that we can find several classes of 
solutions from a totally algorithmic procedure without using additional constraints but 
proceeding only with the considered PDE system. We made an analysis of the symmetries 
of the system ([1]) modelling a planar flow of an ideal plastic material and we computed 
its infinitesimal symmetry generators. Subsequently, using these symmetries, we applied 
the SRM to find G-invariant solutions of this system and this led to several solution 
classes. Many types of solutions were found. For example we got rational algebraic fl37p . 
trigonometric (157|) . inverse trigonometric (1191) . implicit (155]) and some solutions in term 
of one or two arbitrary functions of one variable (see ( I114p and (I122p ). which can be 
chosen to be Jacobi elliptic functions. Some others obtained solutions are expressed in 
term of quadratures (l52l) that must be solved numerically. However, it is more simple to 
numerically solve these quadratures than to completely integrate the reduced equations 
numerically. 

The application of the method discussed above yielded several families of tools that 
can be used for the extrusion of ideal plastic materials. For each solution family one 
may draw extrusion dies corresponding to compatible choices (with the flow lines of the 
solution for the velocities u and v) for given angles and feeding velocities. Since we are 
free to choose the parameters in the solutions and to select which lines of flow to use, 
we consequently have a wide variety of tools for each class of solutions. One should note 
that the similarity solutions corresponding to the additive separation of the velocities 
is particularly interesting since it is expressed in term of arbitrary functions and this 
enlarges the variety of possible tools. In this paper, we have drawn some examples of 
possible extrusion dies. For example we traced a tool that could curve a material in the 
figure [Hand one in the shape of a deformed cross admitting two mouths to feed the tool. 
A tool that could be used to thin a plate is shown in figure [3] and one to wave a plate 
in figure HI A particularly interesting tool from the applications standpoint is shown in 
figure O It might be used to shape a ring by extrusion. 

The method of characteristic^ has already been used to obtain some solutions of the 
system ([1]). However, with this method one is constrained to require the curves limiting 
the plasticity region to be characteristic curves. This is not, in general, the case for G- 
invariant solutions. Generally, when we use the method of characteristics to find solutions 
of ([I]), we obtain u and v by numerical integration along the characteristics (as in-), while 
the use of the SRM ables us to find some classes of analytical solutions not covered by 
the method of characteristics. In this context, a natural question arise: what physical 



34 



insight does one gain from exact analytic particular solutions. A partial answer is that 
they show up qualitative features that might be very difficult to detect numerically: the 
existence of different types of periodic solutions and different type of locahzed solutions. 
Stable solutions should be observable and should also provide a good starting point for 
perturbative calculations. 

The next step of this work is a systematic Lie group analysis of the symmetries of 
the nonstationary system in a 2 + 1 dimensional modelling plane flow of an ideal plastic 
material. We expect that the SRM will lead to wider classes of physically important 
solutions and consequently to new extrusion dies. 
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Table II: List of one-dimensional subalgebras and their invariants. To ease the notation, we noted the invariants Fi = F, F2 = G, F3 = H. 



The index i in Cij corresponds to the subalgebra dimension and the index j to the number of the subalgebra. 









symmetry 






Invariants 






No 


subalgebra 


parameters 


variable ^ 


F 




G 


H 




£1,1 


{aPi + a«(l - <5oa)P2 + flfa} 


a ^ —a 


<I™(<5oa ~ + "^2/ 


u 




V 


a — 


a+a"(l— 5oa) 




f 1 if a = 

where <5oa = 

[0 if a ^ 


a = —a 


2/ 


u 




V 


(T — 


fix 






„-)P,i a 7^ -a. 

a = 1 = —a 


a™a; — ay 
x + y 


u a ''+\ 

—X + U 




y {l-oi)(x+y) 
OL-\-a.^ 

V 


a — 
a — 


" CK + O 

X 






a / 1 / -a 


a™x — ay 




— 2oixy 


^, a(x^-y-')+2a°'xy 




P{x+y) 


r t > 




U 1 -i 




0, -I- a 


a — 








a = 1 = — (/ 


,r + y 


-r ( :^ T y ) + 




^ + '; 


<j — 


„J,r 


-^1,4 


{D2 + aPi +a"P2 + 6P3} 


a ^ —a 


— a"a; + ay 

x + y 






ve 


c — 
c — 


ba; 


■^1,6 


{Di+a«D2+6P3} 




y. 

X 








(T — 


bin a; 


•Cl,6 


{Di + D2 H-B + bPa} 




y 

X 


X X 




^-Inx 

X 


(T — 


bin a; 




+ aPa + aP4 + b^Ps} 




u. 

X 


u — a\nx 




V — b°'lnx 


(T — 


a In a; 


•Cl,8 


{aP4 + a«P5+/3P3} 




?i = a:, 6 = y 


(a + a«)CT - /3 


u + v) 


av — a"u 






-^1,9 


{£12 + bPa} 




ii=x, i2=y 


(T — b In u 




u 

V 






-^1,10 


{B + aPa} 




^i=x,^2=y 






-u + v 

V 







Table III: List of two-dimensional subalgebras and their invariants. To ease the notation, we denoted the invariants Fi = F, F2 = G, 
F3 = H. The index i in Cij corresponds to the subalgebra dimension and the index j to the number of the subalgebra. 









symmetry 




Invariants 








No 


subalgebra 


parameters 


variable ^ 


F 


G 




H 




■^2,1 


{Pi + aP2 + aP4, P3I 




ax + y 




—ax + u 


V 












a + a" 7^ 


-a'^x^ 


- ay 


-b'^u + uiv 




d {u-\-v) 










a ^ I ^ —a, 
a; = -b = 1 




- <^y 


U + V 


d(x + v) 


2~ 










a = —a = 1, 
a; + b"^ 7^ 


— -1 




-b'^u + UJV 




u-\-v 
uj + b" 










a = —a = 1, 
o) = -b= 1 


-a"a;-f 


- cty 




" P 2 


u — v 
~2~ 






^2,3 


{B + aPi+ a"P2, /3P4 + bfP5 + P3} 


a" 7^ 




- ay 


x^ 

a. 1 

2 

„2 

- (1 - ajxy + — (/3 + bi^ja 

Z 


tt 1) 1 r\i 


2 ^ ' y 

- b^)(j 










a = 1, 
a = 


y 




u + v- ^ +xy + {p + bl^)a 




-xy + {0- b^)a 






£■2,4 


{D2 + aPi + a"P2 + bPa, P5} 


a + a« 7^ 


-a"x-{ 


- ay 


x+y 

ue o+o™ 






c — 


b{x+y) 
a+a^ 


a = —a = 1 


x + y 




ue 2 


x — y 








'C-2,5 


{aPi + a«P2 + pPz,ojPi + b-^Ps} 


MO, 

a + a" 7^ 


-a^x -\ 


- ay 


we o+n" 






(J — 


b{x + y) 
Q + a'^ 






MO, 


x + y 




lie c»+£i'^ 






(T — 


b{x-y) 






a = —a = 1 














2 


C-2,& 


{B + aPi + o«P2, /3P4 + b'^Ps} 


b7^ 




- ay 


2 - b'^v^ + a (^(b'^a - I3)x^ 
+ (1 - q) (^-b"y^ - 2(3xy) 


- 2b'^xy'^ 




a 




C-2,1 


{B + aPi + a"P2 + Pa.m + b'^Ps} 


b^O, 
a« 7^ 


—a" 3; - 


- ay 


I3v - bl^u - al31^ - bf' ^ - (1 


- a)l3xy 




a — 


V 






MO, 

a = 1, a = 


y 




- b^u -13^- bl^xy 






u — 


X 



Table IIIII continued 



•^2,8 


{B + aPi + 


a°'P2 + P'MliPi + hf^Ph} , 


fe ^ 0, 
a + a" 


7^0 


—a^x + ay 




V — a^ fl — a)xv 


a 


x + y 








b^U, 
a = l, 


a = -1 


X + y 




x^ 

2 


a 


2 


-^2,9 


{B + aPi + 


a°'P2 + P3,l3P4 + bf^P5} 






b^li 

— CL^ X + ay 




V — a^ (1 ^ a)xy 


a 


x^-y 


-^2,10 


{B + aPi + 


a°'P2 + P3,PP4 + b^P5} , 


b^O 
a + a" 


^ 


—a^x + ay 




V 


a 


x-\-y 








b^O, 
a = l, 


a = -1 


X + y 




V 


a 


2 


^2,11 


{D2 + aPi - 


^a°'P2+bP3,Pi + cP5} 


a + a" 




—a°'x + ay 




x + H 

(v — c«)e °+"" 


a 


- h^Y^ 






a = 1, 


-a = 1 


X + y 




(v — CM)e 


a 


b ^ 

OL — a 


-^2,12 


{Di +D2 + 


bP3,B + P3\ 






y 

X 


cr + - - b\ny 


a - ^ - blnx 

X 






-^2,13 


{Di +D2 + 


bPi,B + P-i\ 






R 


u — a\nx 


V — \rvx 






C2,14 


{D1+D2 + 


bP3,B + P3\ 









u — \i\x 


hv -\- ahix — a 






■^2,15 


{Di + D2 + 


bPs,B + Ps} 









u — £a + ealiix 


V — bu — \nx 






-^2,16 


{Di + aD2 


f feP3,P3} 









ux^"- 


vx~^ 






C2,17 


{DiaPz,D2 


+ ^^3} 








-r 


blnu + alnx — a 


V 

11. 






-^2,18 


{Di +D2 + 


B,P3} 






y 


u + i/ In a; 


v-\-x\i\x 














X 






■^2,19 


{Di+aPz,D2\ 








-r 


—?r^ 

11. 




a 


— a In X 


•^2,20 


{Di + aD2 


^bP3,Pi + cP5} 


c^O 




y 


(cu — v)x~'^ 




a 


— blnx 


^2,21 


{Di +D2 + 


B + aP3,Pi + bP5} 


b^Q 




y 


bu — V ''r {x by ) In x 




a 


— a\nx 


^2,22 


{Di + aD2 


i-bP3,B} 






y 

X 






a 


— blnx 


^2,23 


{Di + aPi - 


i-a^-Ps+bPi^B} 






y 


— {ax-\-a^' y) In x-\-xu-\-yv 






— blnx 






X 


y 




a 


^2,24 


{Di+P4+aP3,P5} 






y 


u — a\nx 




a 


— blnx 


-^2,25 


{Di+P5,P4+aP5} 






S 


V — bu — In X 




a 


— a\nx 


•^2,26 


{Di + aD2 


f feP4,P4} 













a 


— blnx 


•C2,27 


{Di + D2 + 


B + aP3,P4} 








X 




^ - Inx 

X 


a 


— a\nx 



Table IV: List of one-dimensionnal subalgebras and their invariants. To ease the notation, we denoted the invariants Fi = F, F2 = G, 
Fs = H. The index i in Cij corresponds to the subalgebra dimension and the index j to the number of the subalgebra. 









symmetry 






Invariants 


No 


subalgebra 


parameters 


variable ^ 


F 




G H 


■^2,28 


{Pi + aP2 + bPa, P2 + P4 + CP5 + aPa} 




{aa — b)x — ay + a 


ax — y + u 




{ax — y)c + V 


-^2,29 


{Pi + aP2 + P3, Pi + bP2 + P4 + cPs} 




^''7 + a 

a—h 


u + -"^+^ 

a—o 




1 c(-ax+y) 
a—h 


■^2,30 


{Pi+aP3,P2 + bPa} 


+ a 7^ 


a — {ax — ay) 


U 




V 


-^2,31 


{D2 + aPi + (1 - a)P2 + aP3, 
(l-a)Pi+6i-"P2 + P3} 




u — a[ax + y) 
-(1-a) 

• ((l-a&i-«)a; + ay) 






""^) yg{-(^-°')i-b^~''x+V)-°'x) 


-^2,32 


{DiaP4 + a"P5 + &P3, ^Pi + c^P2} 




a — a\Yi{xb^ — (3y) 


u—x 




v-f3bx+{P-l)y 






xb^ —0V 


■^2,33 


\Di+D2 + bP3,l3Pi+cl'P2\ 




a — 61n(a;c^ — Py) 


u — a ln(c"a; — 




V — a" ln(c^x — Py) 


■^2,34 


|P>i + aD2 + bPz, /3Pi + 6'^P2} 


67^0 


a — bln{c^x — jiy) 


u{d^x- iSyy 




v{cf^x — Py)~'^ 



Table V: Subalgebra sort by the form of the invariant solution. The functions /, g, <j), 7 depend on x, y and F, G depend on ^. 







Ansatzs 


No. 




u = } + t\ 
v = g + G 


u = fh\ 
v = gG 


u = f + <pb; 
V = g + jG 


u = h{y) + /F 


u = Ghiy) + /F 


u = h{y) + f + (l>F 


V = q{u) + g + 7G 


V = gq{u) + G 


1 


a\x + a2y 


l^l,l-,l^l,2, 
Cl,Z,C2,l 
^2,2, £2,3 






•^2,4, C,2,i 


■^2,6, ^2,7 




-^2,8, 'C2,9 
•^^2, 10, ^2,11 




2 


y 

X 


1^1,7, 1^2,13, 
C-2,U, C2,15 


^1,5, ^2,16 
^2,17 


^1,6, ^2,12, 

^2,18 


C-2,20 


/^2, 22, £2,23 


l^\,%,C,2,2\,C,2,2A 


■^^2,26, ■i^2,26 
■^2,27 


-Ci.io 


3 


aix + a2y + cr 


l--2,2S, C2,29 
/^2.3() 


£2, 30 














4 


as ln(aia; + a2y) + a 


^2,33 


^2,34 


i^2,32 
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